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LISTING OF THE CLAIMS 

1 . (Currently Amended) A computer implemented method for predicting the structure of a 
membrane-bound protein having a plurality of a helical regions, the method comprising: 

identifying a range ranges of amino acids in an amino acid sequence of the membrane- 
bound protein as transmembrane regions of the membrane-bound protein; 
constructing each of two or more helices in a set of helices for the transmembrane 
regions; 

obtaining an optimized structure for each of the two or more helices; 
assembling the optimized structures of the two or more helices into a helix bundle: 
optimizing a helix bundle configuration for of the s e t of helic e s helix bundle using a first 
molecular dynamics simulation; 

after optimizing the helix bundle configuration, constructing one or more inter-helical 
loops to generate a full-atom model of the membrane-bound protein; and 
optimizing the full-atom model using a second molecular dynamics simulation, th e r e by 
providing to provide a predicted structure for the membrane-bound protein. 

2. (Cancelled). 

3 . (Currently Amended) The method of claim 1 , wherein the constructing each of two or 
more helices in the set of helices for the transmembrane regions includ e s one or more of: 
constructing are each of two or mor o canonical right-handed a-h elices corresponding to the 





transmembran e regions, calculating 
helices, and optimizing each of the c 




energy configuration for each of the canonical 



4 - 34. (Cancelled) 
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35. (Currently Amended) The method of claim 1 , wherein the optimizing a helix bundle 
configuration of the helix bundle includes on e or mor e of : 

assembling a h e lix bundle including e ach of th e set of helices, and 

calculating a minimum-energy configuration for the helix bundle in a lipid bilayer. 

36. (Previously Presented) The method of claim 1, wherein: 
the membrane-bound protein is a G-protein coupled receptor. 

37. (Currently Amended) The method of claim 1 , wherein: 

identifying a range ranges of amino acids in the amino acid sequence as transmembrane 
regions includes aligning the amino acid sequence with an experimental or theoretical 
helical template. 

38. (Currently Amended) The method of claim 1 , wherein: wherein identifying a range of 
amino acids in the amino acid sequence as transmembrane regions includes includes: 

determining the periodicity of hydrophobic residues in the amino acid sequence; and 
optimizing a helix bundl e configuration includes identifying a plurality of lipid-accessible 
residues based at least in part on the determined periodicity. 

39. (Currently Amended) The method of claim 1, wherein: 

c on s tructing obtaining an optimized structure for each of the two or more helices in a set 
of helices for the transmembrane regions includes optimizing each of the two or more 
helices in th e set of helices using a torsional molecular dynamics method. 

40. (Previously Presented) The method of claim 39, wherein: 

the torsional molecular dynamics method uses the Newton-Euler Inverse Mass Operator. 



41. 



(Currently Amended) The method of claim 1 , wherein: 
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constructing each of two or more helices4n a set of h e lices for th e transm e mbran e regions 
obtaining an optimized structure for each of the two or more helices includes determining 
3-D coordinates that defin e the structure of for each of the two or more h e lix in the s e t of 
helices. 

42. (Currently Amended) The method of claim 1 , wherein: 

optimizing assembling the optimized structures of the two or more helices into a helix 
bundle configuration includes determining a rotation and tilt of each helix in the set of 
helices. 

43. (Currently Amended) The method of claim 1 , wherein: 

optimizing assembling the optimized structures of the two or more helices into a helix 
bundle configuration includes orienting th e helix axes of the two or more helices 
according to the 7.5 A electron density map for rhodopsin. 

44. (Currently Amended) The method of claim 38, wherein: 

optimizing assembling the optimized structures of the two or more helices into a helix 
bundle configuration includes orienting the identified lipid-accessible residues to face the 
outside of the helix bundle. 

45. (Previously Presented) The method of claim 1 , wherein: 

the first molecular dynamics simulation is a rigid body molecular dynamics simulation. 

46. (Currently Amended) The method of claim 1 , wherein: 

optimizing assembling the optimized structures of the two or more helices into a helix 
bundle configuration for the set of helices includes modeling an effect of an environment 
of the membrane-bound protein , wherein the effect of the e nvironment is simulated with 
a continuum description of a water environment and a lipid bilayer. 
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47. (Currently Amended) The method of claim 45, wherein: 

the first molecular dynamics simulation uses [[the]] a DREIDING force field, charges 
derived from charge equilibriation to simulate lipids in the membrane, and charges from 
CHARMM22 for the membrane-bound protein. 

48. (Currently Amended) The method of claim [[ 1 ,]] 47, wherein: 

the second molecular dynamics simulation is a mixed mode molecular dynamics 
simulation. 

49. (Currently Amended) The method of claim 48, wherein the mixed mode molecular 
dynamics includes : 

the s e cond mol e cular dynamics simulation uses a torsional molecular dynamics m e thod 
to model modeling the helices and inter-helical loops with a torsional molecular 
dynamics method; and 

a rigid body mol e cular dynamics method to model treating lipids in the [[a]] membrane in 
which th e membrane bound protein is situat e d as rigid bodies, and counterions Na + and 
CI" as free cartesian atoms: 

simulating the outside of the lipids with surface-generalized Born model continuum 
solvent description: 

performing constant temperature dynamics with Hoover algorithm for 50 ps with time 
steps of 1 and 5 fs: and 

using a dielectric constant of 60.0 to simulate the low dielectric region surrounding the 
membrane . 

50. (Currently Amended) The method of claim 1 , wh e r e in: wherein the second molecular 
dynamics simulation includes; 

dynamic optimization of the structure using cell multipole methods for calculation of 
nonbond forces, and [[or]] 
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fast torsional dynamic methods selected from Newton-Euler Inverse Mass Operator and 
Hierarchical Newton-Euler Inverse Mass Operator . 

5 1 . (Previously Presented) The method of claim 1 , wherein: 

at least the second molecular dynamics simulation includes a solvent approximation. 

52. (Previously Presented) The method of claim 51, wherein: 
the solvent approximation is a continuum solvation model. 

53. (Previously Presented) The method of claim 52, wherein: 

the solvent approximation includes the Surface Generalized Born model or the Poisson- 
Boltzmarm description. 

54. (Previously Presented) The method of claim 53, wherein: 

the solvent approximation is an empirical approximation comprising estimating solvation 
free energy as a function of solvent accessible protein surface area. 

55. (Previously presented) The method of claim 1, wherein: 

the second molecular dynamics simulation is performed for a time in the range from 
about 1 00 ps to about 1 ns. 

56. (Previously presented) The method of claim 1 , wherein: 

the set of helices includes four or more membrane-spanning a-helices. 

57. (Previously presented) The method of claim 1 , wherein: 

the set of helices includes seven membrane-spanning a-helices. 



58. 



(Cancelled) 
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59. (Currently Amended) The method of claim [[50,]] L wherein the optimizing the full- 
atom model further includes : 

prior to the second molecular dynamics simulation, 

performing a fall atom minimization of the fall-atom model with a barrel of lipid 
surrounding the protein is p e rformed . 

60. (Previously presented) The method of claim 1 , wherein the amino acid sequence of the 
membrane-bound protein is obtained from GeneBank. 

6 1 . (Previously presented) The method of claim 1 , wherein the predicted structure is output 
in protein data bank format. 

62. (Previously presented) A programmable digital computer, configured to perform the 
method of claim 1 . 

63. (New) A computer program product tangibly embodied in a machine-readable storage 
device, wherein the computer program includes instructions for executing the method of claim 1 
on a programmable processor. 



64. (New) The method of claim 45, wherein the rigid body dynamics is carried out for 1 50 
ps. 



